Organized Current Patterns in Disordered Conductors 
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We present a general theory of current deviations in straight current carrying wires with random 
imperfections, which quantitatively explains the recent observations of organized patterns of mag- 
netic field corrugations above micron-scale evaporated wires. These patterns originate from the most 
efficient electron scattering by Fourier components of the wire imperfections with wavefronts along 
the ±45° direction. We show that long range effects of surface or bulk corrugations are suppressed 
for narrow wires or wires having an electrically anisotropic resistivity. 
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Electron scattering by microscopic structural imper- 
fections in thin conducting films is a major factor deter- 
mining their conductivity properties, especially at low 
temperatures^. Orginary polycrystalline metal wires 
with straight boundaries are usually considered to have 
ohmic conductance with a homogeneous current flow on 
a scale much larger than their grain size (typically tens of 
nanometers) 3 . Measurements using ultracold atoms as a 
highly sensitive probe to minute changes in the magnetic 
field have revealed directional deviations of the current 
flow far from the edges of the wire 4 . Recent observa- 
tions of atomic density fluctuations a few microns above 
wires of different thickness and grain size have revealed 
organized patterns of current flow directional deviations 
which are oriented predominantly at ±45° relative to the 
wire axis, on a length scale as large as tens of microns 5 . 
It was shown that this effect is a general property of 
electron scattering by random imperfections in the con- 
ductor. In contrast to previous observations of atomic 
density fluctuations above current carrying wires&£&^£, 
which were attributed to current irregularities due to wire 
edge corrugation o 10 ! 11 ! 12 , the recent observations empha- 
size the importance of the wire surface or bulk structural 
imperfections on a length scale of the order of a micron 
or longer. 

Here, we present a detailed model for the current ir- 
regularities formed in a current carrying wire with ran- 
dom geometrical perturbations or bulk resistivity inho- 
mogeneities. This model enables not only the quanti- 
tative understanding of the observed patterns and their 
origin, but also provides predictions of electron trans- 
port properties in wires with various geometries and crys- 
talline structures. Together with further measurements 
using the ultracold atomic probe, it is expected to shed 
new light on electron transport, and to allow for less cor- 
rugated atomic traps and guides to be developed for atom 
optics and quantum technology^. 

The results of our calculations are demonstrated in 
Fig. [TJ comparing the measured atomic density patterns 
with patterns calculated by assuming random imperfec- 
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FIG. 1: Measured^ and simulated patterns of atomic density 
(normalized) ^3.5/xm above (a-b) a thick wire (H = 2/xm), 
and (c-d) a thin wire (H = 280nm), with grain size of ~ 70nm 
and ~ 40nm, respectively. Atomic density variations measure 
corrugations in the magnetic field component SB X along the 
wire axis x . These corrugations are due to current directional 
deviations from the main current along x. The trapping po- 
tential ensures a cloud width of ~ 1/xra along y and z, and 
hundreds of jim in the x direction. To create the above 2d 
maps, the different transverse locations are scanned across the 
wire at a constant height. The simulation assumes a combi- 
nation of bulk resistivity inhomogeneity and geometrical per- 
turbations of the wire with parameters chosen such that the 
power spectrum of the magnetic field along x, averaged over 
y and over many realizations of the simulation, fits the mea- 
sured power spectrum (right). 



tions of the wire geometry or bulk resistivity. Each spec- 
tral (Fourier) component of these imprefections is a plane 
wave with a random phase and an amplitude taken from 
a non- white isotropic power spectrum modeled in Ref. 0. 
The angular preference of the patterns emerges from a 
universal electron scattering mechanism described below, 
where the apparent difference in the spectral composi- 
tion of the density fluctuations above the thick wire of 
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FIG. 2: (a) A single Fourier component (plane wave) of the 
planar resistivity perturbations due to bulk inhomogeneity 
or wire thickness variations induces current flow directional 
changes (arrows). The current tilts along low resistivity wave- 
fronts and across high resistivity wavefronts. (b) The ampli- 
tudes of transverse current component SJ y (n) [k = (k x ,k y )], 
generated by resistivity perturbations Sp(n) oc 1/k (see 
Ref. 5), are proportional to sin 20k (color scheme). For a 
wire with a finite length L and width W, these amplitudes are 
calculated at discrete values of which are integer products 
of 2tv/L and 2tv/W. (c) As a consequence of this discrete- 
ness, the corrugations at long wavelengths (k x W < 1) are 
suppressed when the wire becomes narrower: the density of 
k y states becomes lower and no corresponding values of k y 
exist along the maximum scattering amplitude line k y ~ k x . 
Here, the power spectrum of the magnetic field corrugations 
at 3.5pm above the center of the wire is shown as a function of 
k x for resistivity perturbations Spn/po = 3.4-10 _4 (^o/^) with 
Ko = 27r/680pm _1 . (d) For comparison, the magnetic field 
corrugations above the center of the wire are shown for wires 
of different widths for a model assuming edge fluctuations 
with \Sy±(k x )\ = 10nmx(Ko/k x )]. Here, short wavelength 
components are suppressed when the wire becomes wider. 

thickness 2 pm (Fig. [T^L-b) and the thin wire of thickness 
280nm (Fig. Hfc-d) is attributed to the different nature of 
the wire imperfections. 

Our calculations are performed for a metallic wire as 
in Fig. [2^l, having a rectangular cross-section of width 
W in the y direction and thickness H in the z direction, 
and carrying a current density J(r) = Jox + SJ(r). Here 
Jo = I/WH is the regular current density for a total 
current /. Utilizing Ohm's law J = E/p where E is the 
electric field, and Maxwell's equation V x E = 0, we find 

V x J = x J , (1) 

P 

with p = po + dp( r ) being the isotropic resistivity. Al- 
though significant resistivity perturbations may exist in a 
polycrystalline metal near grain boundaries (length scale 
of nanometers), one may safely assume that over most 
of the length scales of interest (microns) Sp <C p, such 



that in the Fourier expansion Sp(r) = J2k P^ e ^ r one nas 
|Pk| <C po for any relevant wave number k = (k x ,k y ,k z ). 
By keeping only terms up to first order in the resistivity 
gradient Vp and using the current continuity equation 
V • SJ = we obtain the solution for the components of 
the current irregularities as a function of the bulk inho- 
mogeneity 

5j( bulk )(k) = J (|j ¥ k-x)^, (2) 

where the transverse components (k y , k z ) of the wave vec- 
tor k take the discrete values 2i:{m/W 1 n/H) with inte- 
gers m and n, — oo < m, n < oo. 

The horizontal transverse current irregularities S Jy bulk ^ 
are proportional to k x k y /k 2 oc sin2#k, where #k = 
tan -1 (k y I 'kx) is the angle in the x-y plane. This im- 
mediately implies that transverse currents are predom- 
inantly generated by Fourier components of the resis- 
tivity perturbations with wavefronts oriented at ±45°. 
Vertical current irregularities SJ z hulk ^ are proportional 
to k x k z /k 2 oc sin 20k, where 0k = tan -1 (k z /k x ). Ver- 
tical currents are therefore significant only for Fourier 
components satisfying k x ~ k z (0k ~ ±45°), namely, for 
longitudinal wavelengths 2tt / k x of the order of the thick- 
ness H or less, corresponding to non-zero values of k z . 
For thin wires, these wavelengths are usually beyond the 
spatial measurement resolution in the x — y plane. At 
wavelengths of interest, much larger than H (k x <C k z \ 
vertical currents are suppressed as SJ^ 1 ^ oc <C 1. 

In the following we refer to the spectral regime 
(k x ,k y ) <C 2tt/H as the "thin film limit", where only 
contributions from Fourier terms with k z = are impor- 
tant. We will then consider the film as two-dimensional 
and characterize it by the real-space vector ^ = (x,y) and 
Fourier space vector k = (k x ,k y ). Thickness variations 
of the wire 5H(x,y) may then be regarded as irregular- 
ities of the thin film resistivity £p thlckness = —p SH/H. 
Figure [2^ demonstrates the generation of periodic hor- 
izontal current directional deviations due to resistivity 
perturbations originating from bulk or thickness varia- 
tions. 

A typical magnetic potential along an elongated trap, 
such as that used in Ref. [5|, is determined mainly by the 
longitudinal component of the magnetic field fluctuations 
at the trapping position. Its Fourier spectrum at a height 
zo is related to the current irregularities in the wire by 

/o 
dz' e-^~ z '\ x 
-H 

x[5J y (k x ,ky,z) +ism9 K ,5J z (k x ,ky,z')} , (3) 

where po is the permeability of the vacuum. Here, 
the free space (continuous) Fourier transformations 
SJ (k x , ky , z') may be approximated by their discrete form 
as in Eq. (0) if kW > 1 and z <C W/2. Substitut- 
ing SJ of Eq. ([2]) into this expression one finds that 
SB x (k) oc e -K2;o sin2#ft limiting the spatial resolution 
in the x-y plane by the measurement distance zq. 
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The sin 20 n dependence together with a ~ 1/n de- 
pendence of the resistivity perturbations (found in the 
spectral analysis of the data in Ref . d) , demonstrated by 
the color map in Fig. [2}d, describe well the behavior in 
the continuum limit W — > oo. However, for finite widths 
(and a finite measurement length L) k x and k y assume 
only discrete values which are integer multiples of 2tt/L 
and 2tt/W, respectively, as demonstrated by the grid of 
dots superposed on the color map. It follows that for 
small values of k xi no counterparts k y exist on the grid 
which lie in the region where | sin 26^ | is large, or more 
specifically, around the line K = 45°. This implies that 
at wavelengths larger than the wire width W the cur- 
rent irregularities are significantly suppressed beyond the 
suppression caused by the reduction of grid points. This 
prediction is demonstrated in the power spectrum shown 
in Fig. [2k. This result is very different from the effect of 
current irregularities due to edge roughness, which char- 
acterized measurements of atomic density fluctuations in 
some previous wor kj 1Q i n i 12 . In that case, the short wave- 
lengths are exponentially suppressed near the center of 
the wire, while only wavelengths of the order of the wire 
width or more are effective (Fig. [2]i). 

Another prediction of our model is obtained when we 
generalize the situation to the case where the conducting 
wire is electrically anisotropic, such that the resistivity is 
a diagonal tensor and Ohm's law generalizes to Ej = pjjj 
for j = x, y, z. In this case Eq. ([2]) becomes^ 



£J (bulk) (k) = Jo 



k x 



x 



Px,0 



(4) 



where q = (k x / p x ,k y / p y ,k z / p z ). In the limit of a thin 
film, where k z = 0, the horizontal transverse current 
irregularities SJ y are proportional to sin 20^/(1 + (r — 
l)cos 2 6^), where r = p y /p x is the resistivity ratio. As 
demonstrated in Fig.[3j the scattering at angles On < 45° 
is suppressed if r > 1 and enhanced if r < 1, thus chang- 
ing the preferred scattering wavefront angle in the range 
0° < On < 90°. The overall magnetic corrugations are 
suppressed as r -3 / 4 in the limit of high anisotropy r ^> 1. 

Now we turn to a more detailed theory of current irreg- 
ularities due to geometrical imperfections of the wirei^. 
We solve Eq. (pQ) with Sp — > and with boundary con- 
ditions ensuring that the current flows parallel to the 
boundaries. Taking the upper and lower surfaces of the 
wire at z = ±H/ 2 + Sz± and the right and left edges at 
y = ±W/2+5y±, where 5z±(x, y) and 5y±(x, z) are small 
fluctuations of the corresponding surfaces, we obtain the 
following boundary conditions, 



W 

8J z (x,y,± — ) 



J, 



d5y± 



J, 



1 dx 
d5z± 
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dx 



(5) 



where terms of second or higher orders in Sz± and Sy± 
were omitted. The current irregularities are then written 
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FIG. 3: For a current flowing through an electrically 
anisotropic wire, the perturbation wavefront angle #J£ ax , giv- 
ing rise to maximum transverse electron scattering, will de- 
pend on the ratio r — p y /p x between the transverse and lon- 
gitudinal resistivities (main plot). The two-dimensional maps 
show the predicted atomic density above a wire similar to that 
presented in Fig. [1] in the extreme cases On =0° (r < 1), 90° 
(r ^> 1) and the isotropic case On = 45° (r = 1). Bottom 
inset - the magnetic corrugation amplitude as a function of 
On, which is suppressed when r > 1. For details see Ref. El 



as a sum of two terms, 

5J ( S u rf)(r) = J2e ik ° x [5JZ(y,z) + 6Jl(y,z)] .(6) 

Eq. (pQ) with V p = together with the continuity equa- 
tion V • J = imply that the current can be written as 
the gradient of a potential function Jj( surf ) = VJF, which 
satisfies the Laplace equation V 2 T = 0. It follows that 
the terms in Eq. (j6]) have the form 



53l(y,z) = ik x J h m , ± (k x )e i2 ™ y/W e-^ H W™,(8) 



where the exponential terms describe the attenuation of 
current fluctuations induced by each boundary pertur- 
bation at a distance A n = [k% + (27rn/i^) 2 ] -1 / 2 from the 
left/right boundaries and A m = [k 2 x + (2nm/W) 2 )- 1/2 
from the top/bottom boundaries. Since S J is derivable 
from a scalar function, it follows that each of the vecto- 
rial coefficients a n? ± and h m ^± can be derived from the 
corresponding scalar coefficients. Linear equations are 
obtained for these scalar coefficients when £J of Eq. (j6]) 
is substituted in the boundary conditions (|SJ^£. 

Next, we describe the solutions of these equations for 
a few typical simple cases. The term (Eq. 0) is 

significant when the edge perturbations 5y± are large and 
the measurement height is comparable to the wire width 
W. This situation was discussed in previous work o 11 ! 12 . 
Here we concentrate on the other limit, where the field is 



measured at a low height and a large distance from the 
edges compared to A n such that SJ^ °° for most values 
of k x . The surface height fluctuations 5z^ then generate 
the following current irregularities (Eq. dSJ) 

k ^ sinh(ftif/2) 2 

A /(surf)/ n TnS^P iK €ik CQsh (^) X.mean /,qX 

^ W - J °Z^ e ^cosh(^/2)^ ' (1Uj 

where SH^ = Jz^ — Sz^ are wire thickness variations 
and Sz^ ean = (Sz^ + 5z^)/2 are height fluctuations of 
the center of the wire. In the thin film limit kB. <C 1, we 
find kz <C 1 for \z\ < H/2 such that SJ y assumes a form 
similar to Eq. ([2]) with Sp^/po —> —SH^/H 

^(r)«-J o ]Te^sin20^. (11) 

In the same limit, the magnitude of the vertical cur- 
rent component 5J Z becomes negligible, since <5ji surf ^ oc 
k x 5z% ean = k x H(5z% ean /H) <C 5z% ean /H. 

Equation ([TT]) shows that the surface roughness 5z± 
causes current irregularities mainly through the thickness 
variations 5H = 8z+ — Sz-. In thin films, long scale sur- 
face height variations which follow bottom (wafer) 
surface variations Sz- , are not expected to cause signif- 
icant current variations. On the other hand, the effect of 
thickness variations is expected to be more pronounced 
in thin films, where SH/H is larger than in thicker wires. 
The power spectrum of the measured magnetic field pat- 
tern above the thin wire (Fig. []J), which was analyzed 
in Ref. 0, could be explained by a model assuming that 
thickness variations exist mainly at short length scales 
(below ~20/im), while 5z+ ~ Sz- at long length scales. 
This may explain the shorter characteristic length scale 
of the features in the thin wire (Fig. relative to the 



thick wire (Fig. QJi), where the measured surface rough- 
ness was not sufficiently large to account for the magnetic 
field fluctuations even if the bottom surface variations 
were assumed to be uncorrelated with the top surface. 
This analysis implies that bulk resistivity perturbations 
with oc 1/k spectrum could play an important role in the 
thick wire, but they are much smaller in the thin wire. 

To conclude, we have presented a detailed model for 
current directional deviations in thin wires with random 
imperfections on a length scale of the order of a micron 
or longer. These deviations may arise either from bulk 
resistivity inhomogeneities or from geometrical pertur- 
bations of the wire, where the significance of each factor 
depends on the wire thickness and fabrication process. 
In both cases, electron scattering is dominant at wave- 
fronts oriented at ±45° relative to the main current axis. 
The model predicts a strong suppression of long wave- 
length current deviations originating from bulk or surface 
corrugations in narrower wires. Electrically anisotropic 
materials are also capable of significantly suppressing 
these deviations. Such analysis opens the road for mate- 
rial engineering to considerably improve atom optics on 
atomchips where currents are used for creating magnetic 
potentials for atom trapping and guiding. Comparison 
of this theory with further cold atom magnetometry or 
other measurements providing high field sensitivities and 
spatial resolution, will enable deeper understanding of 
electron transport in thin films. 
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tuations in a wire with rectangular cross section. The case 
of top surface fluctuations was given in Ref. [H only for 
symmetric fluctuations that contribute to the magnetic 



field above the middle of the wire. 

Y. Japha, 0. Entin-Wohlman and R. Folman, unpublished. 



